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Preface 


This  report  is  a  collection  of  material  that  has  been  used 
in  courses  on  search,  detection  and  localization  modeling.   Its 
organization  follows  to  some  extent  material  by  S.  M.  Pollock  in 
Selected  Methods  and  Models  in  Military  Operations  Research  which 
is  listed  in  the  report  bibliography.   The  report  is  not  intended 
to  be  a  text  on  these  subjects.   In  particular,  in  some  areas  it 
does  not  provide  the  depth  of  coverage  that  is  found  in  the  book 
Search  and  Detection  by  Alan  R.  Washburn  which  is  cited  in  this 
report  as  Reference  21. 
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I.   Detection  Models  and  Signal  Detection  Theory 

Signal  detection  theory  is  the  basis  for  analyzing  the 
detection  models  that  are  described  in  this  report.   In  signal 
detection  theory,  the  decision  making  portion  of  a  detection 
system  is  called  the  receiver  and  a  detection  experiment  is  the 
observation  by  a  receiver  of  input  data  generated  during  some 
time  interval.   The  data  that  is  related  to  a  target  is  called 
signal.   The  data  that  is  not  related  to  the  target  is  called 
noise.   In  general,  the  target  data  is  associated  with  a 
localization  region  that  in  some  cases  is  called  a  resolution 
cell.   When  a  detection  experiment  is  performed,  either  the  event 
H0  =  {the  receiver  input  is  noise}   or  its  complement 
H^  =  {the  receiver  input  is  signal  and  noise)   will  occur.   In 
the  first  detection  models  that  are  described  here,  after 
analysis  of  the  input  data  by  a  receiver,  either  the  event 
Dq  ={the  receiver  decides  the  input  is  noise}   or  its  complement 
D-^  =  {the  receiver  decides  the  input  is  signal  and  noise}   will 
occur.   Detection  models  for  which  D±      is  the  complement  of   D0 
are  called  binary  detection  models  or  forced  choice  detection 
models.   Four  events  which  are  important  in  binary  detection 
model  are  indicated  in  the  Venn  diagram  of  Figure  1. 

The  Venn  diagram  emphasizes  a  decision  problem  that  is 
associated  with  a  receiver  that  can  be  modeled  using  a  binary 
detection  model.   The  problem  is  this:   Under  what  conditions 
should  the  event   D]^   occur?   That  is,  under  what  conditions 
should  a  receiver  decide  that  the  input  data  indicates  a  target 
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was  present  in  an  observed  localization  region  either  at  the  time 
or  prior  to  the  time  of  the  observation? 


D, 


Hi 


H- 


D0   n   H0 

D0    n    H^ 

DX   n  H0 

Dx  n  Hx 

Figure  1.   Four  events  of  importance  in  binary  detection  models. 


In  the  detection  model  descriptions  that  follow,  the 
following  notation  and  terminology  are  used:   pf  =  PfD^jHo),   the 
probability  of  D±   given  H0,  is  called  the  false  alarm 
probability;   p^  =  PfDjjH!),  the  probability  of   D]^   given  Hlf 
is  called  the  detection  probability  and   P  =  P^-^)  ,   the 
probability  of   H]_,  is  called  the  prior  probability.   It  is  the 
probability  that  a  target  will  be  in  the  localization  region  when 
the  input  data  is  generated. 

In  the  detection  models,  the  input  to  a  receiver  consists  of 
a  sequence   yi»«'«ym  that  is  a  realization  of  a  discrete 
parameter  stochastic  process  whose  random  variables  are  elements 
of  a  continuous  parameter  stochastic  process.   Although  the 
receiver  input  in  some  cases  can  be  considered  to  be  described  by 
a  continuous  stochastic  process,  because  of  the  finite  amount  of 
information  (unique  data)  contained  in  a  bounded  sequence  of 
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a  continuous  stochastic  process,  because  of  the  finite  amount  of 
information  (unique  data)  contained  in  a  bounded  sequence  of 
finite  length,  a  discrete  stochastic  process  is  sufficient  to 
represent  the  input  in  these  cases.   This  is  established  formally 
by  the  stochastic  sampling  theorem.   In  the  models,  for  an 
observation,  the  input  stochastic  process  has  the  following 
characteristics:   It  is  a  noise  process  when  there  is  no  target 
data  and  it  is  a  noise  process  plus  a  signal  process  when  there 
is  target  data.   In  the  first  detection  model  that  is  described 
in  Section  III,  the  signal  process  is  a  deterministic  process. 
In  the  second  and  third  detection  models  that  are  described 
there,  the  signal  process  is  a  random  process.   To  specify  the 
noise  process  or  a  random  signal  process,  one  needs  only  to 
specify  the  joint  distribution  of  the  finite  sequence  of  random 
variables  that  determine  the  process.   If  the  signal  process  is  a 
deterministic  process,  the  signal  values  are  a  sequence  of  values 
that  are  assumed  to  be  known  before  an  observation  is  performed. 
Consequently,  in  this  case,  to  specify  the  process,  one  needs 
only  to  specify  these  values. 


II.   Decision  Criteria 

To  simplify  the  discussion  of  decision  criteria  and 
decision  rules,  a  receiver's  input  will  be  assumed  to  be 
determined  by  a  single  decision  random  variable  Y.   In  this 
case,  the  input  process  in  determined  by  the  conditional 
distribution  function  Fy(y|H0)   when  the  input  is  noise  alone 
and  by  the  conditional  distribution  function  FyCylH]^)   when  the 
input  is  signal  plus  noise. 

The  condition  that  a  receiver's  input  is  required  to 
satisfy  in  order  that  the  event  D±     will  occur  can  be  specified 
in  terms  of  a  decision  rule.   For  the  assumed  case,  a  decision 
rule  is  a  rule  which  determines  for  every  observable  value  of  Y 
the  decision  that  the  receiver  is  to  make.   The  decision  rule  can 
be  considered  to  be  a  function   0(y)   which  relates  each 
observable  value  y  to  one  or  the  other  of  the  following  two 
commands:   d0  =  "Decide  that  the  receiver  input  was  noise."  and 
d1  =  "Decide  the  input  was  signal  and  noise.".   Choosing  a 
decision  rule  0(y)   defines  a  set  ft  of  observable  values  of  Y 
such  that  the  event  D±   =  {  Y  e  ft  } . 

The  problem  which  was  considered  in  Section  I  can  now  be 
restated  in  the  following  way:   What  criterion  should  be  adopted 
in  order  to  determine  a  decision  rule  or,  equivalently,  its 
corresponding  set  ft  ?  A  desirable  characteristic  for  a  criterion 
is  suggested  by  the  following  argument:   Consider  the  odds  in 
favor  of  Hi     given  y   is  observed.   That  is,  consider 


P(H1|Y  =  y)/P(H0|Y  =  y) .   One  might  expect  that  y  would  be  a 
member  of  the  set  n   if  and  only  if  y  made  this  ratio  equal  to 
or  greater  than  some  value  k.   But  this  is  equivalent  to  defining 
ft  as  follows:   ft  =  {  y  :  L(y)  >  K  }  where  L(y)   is  the 
likelihood  ratio  associated  with  an  observed  value  y  and  K  is 
a  constant  related  to  the  constant   k.   This  suggests  that 
choosing  an  optimum  criterion  is  equivalent  to  choosing  an 
optimum  value  for  K.   Four  specific  decision  criteria  are  defined 
next  in  terms  of  K.   For  each  criterion,   ft   has  the  above 
form.   But  for  each  criterion  the  choice  of   K   is  different. 
The  decision  criteria  are: 

1.  The  Neyman-Pearson  Criterion:   Choose   ft   so  that  p^ 

is  a  maximum  subject  to  the  constraint  that  pf  <  a  where  a   is 
a  specified  value.   For  a  continuous  decision  random  variable, 
the  constant  K   is  chosen  so  that  pf  =  a. 

2.  The  Bayes  Criterion:   Choose  ft  so  that  the  expected 
cost  of  a  receiver's  decision  is  a  minimum.   For  a  continuous 
decision  random  variable,   K  =  [ (c10-c00)/ (c01"cll) ] (1-p)/p  if 
c10  >  c00   anc*  coi  >  cll  where  c^j   is  the  cost  of  Dj^  n  Hj . 

3.  The  Ideal  Observer  Criterion:  Choose  ft  so  that  the 
probability  that  the  receiver  makes  an  incorrect  decision  is  a 
minimum.   K  =  (1-P)/P  for  a  continuous  decision  random  variable. 

4.  The  Minimax  Criterion:   Choose   ft  when  P  is  unknown  so 
that  the  maximum  expected  cost  of  a  receiver's  decision  is  a 
minimum.   If  c10  >   c00   and  c01  >  cll'  tnen 


K  =  [ (Cio~cOo)/ (c01"cll) 3 (1-P*)/P*   for  a  continuous  decision 
random  variable.   Here,   P*  is  the  value  of  the  prior  probability 
P  which  would  make  the  expected  cost  of  a  receiver's  decision  a 
minimum  for  this  choice  of  K. 

If  a  model  which  specifies  the  conditional 
distributions   Fy(y|H0)   and  FyCy!^)   and  a  decision  rule  are 
adopted,  then  the  value  of  pf  and  the  value  of  p^  are 
determined.   This  pair  of  values   (Pf,Pd)   i-s  called  a  receiver 
operating  point.   If  the  decision  rule  results  from  using  a 
likelihood  ratio  criterion  such  as  one  of  the  four  listed  above, 
then  it  will  involve  the  parameter  K  through  the  relation 
n  =  {y:  L(y)  >  K} .   And,  for  a  given  value  of  K,   since  n 
uniquely  determines  the  pair   (Pf /Pd) /  a  single  operating  point 
results.   By  varying  K,  a  set  of  operating  points  can  be 
generated  which  determines  a  receiver  operating  characteristic 
curve  or  ROC  curve.   Different  ROC  curves  can  be  produced  by 
changing  either  one  or  both  of  the  conditional  distributions 
which  implies  changing  either  the  signal  process  or  the  noise 
process. 

A  decision  rule  which  results  from  using  a  likelihood  ratio 
criterion  in  a  model  in  which  the  input  process  is  determined  by 
a  set  of  m  random  variables  can  be  expressed  in  terms  of  a  set 
n  as  follows:  n  =  {  (yi,...,Vjn)  •  L(Yi/  •  •  •  /Ym)  -  K  )  where  K 
is  specified  in  the  same  way  as  in  the  case  in  which  m  =  1. 


III.   Three  Binary  Detection  Models 

Three  detection  models  are  examined  in  this  section. 
For  the  first  two  detection  models,  the  input  stochastic  process 
for  an  observation  is  defined  by  a  time  sequence  of  continuous 
random  variables.   The  random  variables  represent  a  sample  from  a 
continuous  parameter  stochastic  process  which  is  sampled  at  times 
such  that  the  random  variables  are  independent.   For  the  third 
detection  model,  the  input  stochastic  process  is  a  counting 
process  and  it  is  defined  by  a  single  discrete  random  variable 
that  is  equal  to  the  number  of  events  that  are  counted  during  the 
observation. 

Model  I:   In  the  first  detection  model,  a  sampled  noise 
value  is  a  value  of  a  normally  distributed  random  variable  with 
mean  zero  and  with  known  variance  a2      and  a  sampled  signal  value 
is  a  known  value  of  a  deterministic  variable.   Thus,  the  input 
process  corresponding  to  an  observation  consists  of  some  number 
m   of  independent  normal  random  variables   Y  ]_,...,  Ym  each  with 
variance   a2.   And,  for  i  =  1,2,3,  •••  ,  m,  when  a  signal  is  not 
present  the  mean  of   Y^   is  zero  and  when  a  signal  is  present  the 
mean  is   s^.   The  result  of  using  a  likelihood  ratio  decision 
rule  in  the  model  can  be  expressed  in  terms  of  a  random  variable 
Z.   This  random  variable  is  called  a  crosscorrelation  statistic 
and  it  is  defined  by   Z  =  Z  s^'Y^   where  the  sum  index   1  =  1,2, 
•••  ,m  here  and  in  the  remainder  of  this  section.   However,  it 
is  more  convenient  to  express  the  result  in  terms  of  a  statistic 
V  which  is  defined  by  V  =  Z/az.   In  terms  of  this  random 


variable,  the  conditional  probabilities  pf   and  p^   are  given 
by:   pf  =  1  -  *(v*)   and  p<j  =  1  -  *  (v*  -  d*)   where   * 
symbolizes  the  standard  normal  cumulative  distribution  function, 

the  threshold  v*  =  (l/az) (a*ln  K  +  (1/2)  2  s[)    and 
d  =  2  sj/a  .   The  parameter  d   is  called  the  detection  index. 
Often,  the  input  stochastic  process  represents  a  quantity 
whose  square  is  proportional  to  power.   In  such  a  case,  the 
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average  receiver  input  power  is  the  random  variable   S  Y-j/m.   If 
a  signal  is  not  present,  the  expected  average  receiver  input 
power  is   N  =  Z  a*/m  =  a*   where  N   is  called  the  noise.   The 
average  receiver  input  signal  power  is   S  =  S  s^/m  where  S   is 
called  the  signal.   In  these  terms,   d  =  m* (S/N)  where   S/N   is 
called  the  signal-to-noise  ratio. 

If  a  receiver's  input  can  be  considered  to  be  a  time 
sequence  of  continuous  voltage  values  such  as  in  the  case  of  a 
sonar  receiver,  in  some  cases  a  frequency  representation  can  be 
used  that  involves  the  concept  of  receiver  bandwidth.    In  these 
cases,  the  noise  process  is  assumed  to  be  such  that  m  =  t/St 
where  t   is  the   integration  time   (the  duration  of  an 
observation)   and   St   is  the   time  between  samples   with 
5t  =  1/[2(BW)]   where   BW   is  the  bandwidth   and   St   is 
determined  by  the  sampling  theorem.   This  implies  that  the 
detection  index  can  be  written  as  d  =  2t» (BW) (S/N) .   By  defining 
N0  as  the  power  spectral  density  where  N0  =  N/BW   the  detection 
index  can  also  be  written  as  d  =  2t« (S/N0) .   In  Reference  2,  the 
conditions  required  by  the  first  model  are  called  Case  I   and  in 


the  following  sections  the  first  model  is  called  the   Case  I 
model.   A  receiver  that  processes  data  such  that  it  would 
implement  a  likelihood  ratio  decision  rule  under  the  conditions 
of  the  first  model  is  called  a  crosscorrelation  detector.    If 
the  description  of  the  input  noise  is  adequate,  a  Case  I  model 
can  be  used  to  obtain  an  estimate  of  an  upper  bound  on  a 
detection  system's  performance,  since  all  the  information 
necessary  to  define  the  signal  is  assumed  to  be  known. 

Model  II:   In  the  second  detection  model,  a  sampled  noise 
value  is  a  value  of  an  independent  normal  random  variable  with 
mean  zero  and  with  variance  a*.      And  a  sampled  signal  value  is 
an  independent  random  variable  with  mean  zero  but  with  variance 
Og.   Thus,  the  input  process  corresponding  to  an  observation 
consists  of  some  number  m  of  independent  normal  random 
variables   Y1,...,Ym  each  with  mean  zero  and  each  with  variance 
a*      when  a  signal  is  not  present  and  with  variance   a*  +   os  when 
a  signal  is  present.   The  result  of  applying  a  likelihood  ratio 
decision  rule  in  this  model  can  be  expressed  in  terms  of  a 
statistic   X  which  is  defined  by  X  =  Z  Y^.   When  a  signal  is 
not  present,  the  statistic   X/N  has  a  chi-square  distribution 
with  m  degrees  of  freedom.   When  a  signal  is  present,  the 
statistic   X/(N+S)   has  a  chi-square  distribution  with  m 
degrees  of  freedom.   So,  in  terms   of  these  statistics,  the 
conditional  probabilities  pf   and  p^   are  defined  as  follows: 
Pf=  p(xm  ^  X*/N)   and  pd  =  P{X^  >  (x*/N) [1/ (1+S/N) ] } 


where   X^   is  a  chi-square  random  variable  with  m  degrees  of 
freedom,  x*  is  a  number  which  is  determined  by  the  decision  rule 
and   S/N   is  the  signal-to-noise  ratio.   A  receiver  that  would 
implement  a  likelihood  ratio  decision  rule  under  the  conditions 
of  the  second  model  is  called  an  energy  detector  (or  square  law 
detector) . 

The  mean  of  a  chi-square  random  variable  with  m  degrees  of 
freedom  is  m  and  the  variance  is   2m.   By  the  central  limit 
theorem,  as  the  number  of  degrees  of  freedom  of  a  chi-square 
random  variable  becomes  large  it  can  be  approximated  by  a  normal 
random  variable  with  the  same  mean  and  variance.   Hence,  for  a 
sufficiently  large  sample  size  m  or  equivalently,  for  a 
sufficiently  large  time  bandwidth  product   t* (BW) ,   the 
conditional  probabilities  pf  and  p^  can  be  approximated  by 
pf  =  1  -  *(v*)   and  pd  =  1  -  ${ [1/(1+S/N) ] (v*  -  d*)}   where  the 
threshold  value  v*  =  (x*  -  mo2)/[(2m)^a2]   and  d  =  t- (BW) (S/N) 2 
if  the  concept  of  bandwidth  is  applicable.   If,  in  addition,  the 
noise   N   is  significantly  larger  than  the  signal   S,  that  is,  if 
t- (BW)  »  1   and   S/N  «  1,   then  pf   and  pd   can  be 
approximated  by:   Pf  =  1  -  *(v*)   and  p<j  =  1  -  *(v*  -  d^). 
In  Reference  2,  the  conditions  required  for  this  approximation 
are  called  Case  II  and  in  the  following  sections  the  limiting 
form  of  the  second  model  is  called  the   Case  II   model.   In 
Reference  2,  an  example  is  given  that  uses  the  Case  II  model  to 
estimate  the  performance  of  a  passive  sonar  system. 
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Model  III:   In  the  third  detection  model,  a  sampled  noise 
value  and  a  sampled  signal  value  are  values  of  independent  random 
variables  that  are  determined  by  independent  Poisson  processes 
observed  for  a  time  interval   t.   The  noise  process  is 
characterized  by  a  counting  rate   a,  the  signal  process  is 
characterized  by  a  counting  rate  as  and  the  noise  and  signal 
processes  are  additive.   Consequently,  when  the  input  is  noise 
alone,  the  input  is  the  value  of  a  random  variable  Y   that  has  a 
Poisson  distribution  with  parameter  at.   And,  when  a  target  is 
not  present,  the  input  is  a  value  of  a  random  variable  that  has  a 
Poisson  distribution  with  parameter   (a  +  as) *t.   For  a 
likelihood  ratio  decision  rule,   pf  =  1  -  P(y*;at)   and 
Pd  =  1  -  P[y*;(a  +  as) «t]   where   P(y;9)   symbolizes  the 
cumulative  Poisson  distribution  function  with  parameter  9   and 
y   is  a  threshold  value  of   y  determined  by  the  decision  rule. 
In  terms  of  the  above  notation,  both  the  mean  and  variance  of  a 
Poisson  random  variable  are  equal  to   9.   When   9   is  large,  the 
cumulative  Poisson  distribution  function  can  be  approximated  by 
the  cumulative  distribution  function  of  a  normal  random  variable 
having  the  same  mean  and  variance.   Using  this  approximation  for 
cases  where   at  is  large  gives   pf  =  1  -  *(v*)   and 
pd  =  1  -  *{[l/(l  +  as/a)*](v*  -  d') }   where  v*  =  (y*  -  at)/(at)* 
and   d  =  a«t«(as/a)2.   In  addition,  if  the  signal  counting  rate 
is  significantly  smaller  than  the  noise  counting  rate,  that  is, 
if   a-t  »  1   and   ag/a  «  1,   then  Pf   and   pd   can  be 
approximated  by:   pf  =  1  -  *(v*)   and  Pa  =  1  -  *(v*  -  d'). 
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This  detection  model  might  be  used  to  describe  a  receiver  whose 
input  for  an  observation  is  the  number  of  photons  counted  by  a 
radiation  detector  in  situations  where   at,   the  expected  number 
of  counts  when  no  signal  is  present,  is  of  the  order  of  thirty  or 
more.   In  this  case,  the  target  might  be  a  wake  segment  for 
example. 

When  a  likelihood  ratio  decision  rule  is  used  in  the  three 
models  discussed  above,  for  the  first  model  and  under  limiting 
conditions  for  the  second  and  third  models,  the  following  result 
is  obtained:   pf  =  1  -  *(v*)   and  p^  =  1  -  *(v*  -  d1)      where  the 
definition  of  v*   depends  on  the  noise  power  N   for  the  first 
and  second  models.   For  a  sonar  receiver  described  by  the  first 
model,  that  is,  by  the  Case  I  model:   d  =  2t« (BW) (S/N) .   For  a 
sonar  receiver  described  under  the  limiting  conditions  for  the 
second  model,  that  is,  by  the  Case  II  model,   d  =  t • (BW) (S/N) * . 
So,  in  either  a  Case  I  model  or  a  Case  II  model  of  a  sonar 
receiver,  the  detection  index  d   is  a  function  of  the  time 
bandwidth  product   t* (BW)   and  signal-to-noise  ratio   S/N.   Since 
sonar  equations  relate   S/N  to  system,  target  and  environmental 
parameters,  a  sonar  equation  can  be  used  to  relate   S/N  to  these 
parameters  in  a  model  of  a  sonar  receiver. 
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IV.  General  Detection  Models 

The  detection  models  that  have  been  considered  to  this 
point  are  based  on  binary  detection  theory.   After  each 
observation,  a  receiver  decides  either  that  the  input 
corresponding  to  the  observation  was  noise  or  else  it  decides  it 
was  signal  plus  noise.   However,  in  some  detection  systems  this 
decision  is  delayed.   In  a  computational  sense,  a  model  of  such  a 
detection  system  is  generally  complex  relative  to  a  binary 
detection  model.   To  illustrate  this,  consider  an  active  sonar 
system  whose  receiver  includes  an  operator.   Suppose  the 
probability  that  the  operator  will  detect  a  target  echo  has  been 
determined  in  a  laboratory  experiment  in  which  the  operator  was 
required  to  decide  after  each  input  corresponding  to  a  resolution 
cell  that  either  the  input  was  a  target  echo  (signal)  and  noise 
or  the  input  was  noise  alone.   In  addition,  suppose  that  under 
operational  conditions  the  operator  normally  delays  this 
decision.   Then,  in  general,  the  probability  that  the  operator 
will  decide  that  the  input  corresponding  to  a  resolution  cell 
that  contains  a  target  is  a  target  echo  and  noise  will  not  be 
equal  to  the  probability  of  the  event  in  the  forced  choice 
experiment.   And,  in  addition,  the  probability  that  the  operator 
will  decide  the  input  corresponding  to  a  resolution  cell  that 
does  not  contain  a  target  is  a  target  echo  and  noise  will  not  be 
equal  to  the  probability  of  this  event  in  the  forced  choice 
experiment.   Consequently,  in  general,  the  value  of  bc'-L   pu   aril 
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Pf   for  an  operational  environment  will  be  different  than  that 
for  the  laboratory  environment. 

One  model  that  has  been  proposed  to  deal  with  this  kind  of 
situation  defines  the  event  that  a  receiver  decides  that  the 
input  corresponding  to  a  resolution  cell  is  signal  and  noise  to 
be  eguivalent  to  the  event  that  out  of   n  consecutive 
observations  at  least   k  of  them  would  result  in  the  decision 
that  the  input  was  signal  and  noise  in  a  forced  choice 
experiment.   The  model  is  said  to  be  based  on  an  k-out-of-n 
detection  criterion.   With  this  criterion,  the  probability  that  a 
target  will  be  first  detected  on  the  jtn  observation  can  be 
found  as  follows:   Determine  the   2^   sequences  of  forced  choice 
responses  that  could  result  for  a  sequence  of   j   consecutive 
observations.   Next,  determine  the  probability  of  occurrence  for 
each  sequence  that  first  satisfies  the   k-out-of-n  detection 
criterion  on  the  j*-*1  observation.   The  probability  of  first 
detection  on  the  jth  observation  is  equal  to  the  sum  of  these 
probabilities.   The  cumulative  probability  of  detection  at  the 
jtn  observation  is  the  sum  of  the  probabilities  of  first 
detection  on  the  i*-*1  observation  for   i  =  l,2,  •••  ,j. 
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V.  Signal-to-Noise  Ratio  Detection  Models 

In  some  radar  and  sonar  detection  models,  for  a 
specified  value  of  pf,  a  minimum  acceptable  value  of  p<j   is 
defined.   This  minimum  acceptable  value  of  p^   and  the  specified 
value  of   pf   define  what  can  be  called  a  minimum  acceptable 
signal-to-noise  ratio   (S/N)m   if   p^   is  a  nondecreasing 
function  of   S/N.   In  some  sonar  detection  models,  this 
signal-to-noise  ratio  in  decibels  is  called  the  detection 
threshold   DT.   In  symbols,   DT  =  10  log (S/N) m.   If  the  minimum 
acceptable  value  of   p^   is  .5,  then   DT   is  usually  called  the 
recognition  differential   RD.   The  difference  between  the 
signal-to-noise  ratio  in  decibels  and  RD   (or  DT)   is  called  the 
signal  excess   SE.   In  symbols,   SE  =  10  log (S/N)  -  RD. 
One  interpretation  of  signal  excess  is  that  for  a 
localization  region  containing  a  target  detection  occurs  with 
probability  one  if  SE  >  0   and  with  probability  zero  if  SE  <  0. 
This  interpretation  provides  the  basis  for  defining  detection  in 
the  three  encounter  detection  models  that  are  discussed  in  the 
Section  VII.   A  more  consistent  interpretation  of  signal  excess 
is:   If   SE  >  0,   then  the  probability  of  detection   p^  is 
greater  than  or  equal  to  the  minimum  acceptable  value   (.5   if 
recognition  differential   RD   is  used  to  define  signal  excess) . 
For  cases  where  p^   increases  rapidly  with  signal  excess  in  the 
neighborhood  of  zero  signal  excess,  the  two  interpretations  are 
essentially  equivalent.   For  a  discussion  of  this  point  as  well 
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as  a  discussion  of  an  operational  case  in  which   receiver 
decisions  are  delayed,  see   Reference  3. 

Signal  excess  (signal-to-noise  ratio)  models  provide  a  basis 
for  general  detection  models,  in  particular,  models  that  describe 
nonstationary  noise  and  signal  processes  and  randomly  changing 
decision  rules.   This  is  illustrated  by  the  models  described  in 
Section  VII.   In  addition,  signal  excess  models  provide  a  basis 
for  delayed  receiver  decision  models.   This  is  illustrated  by  the 
active  sonar  detection  models  in  both  Reference  4  and  Reference  5 
that  are  based  on  a   k-out-of-n   detection  criterion.   In  all  of 
these  models,  the  signal-to-noise  ratio  and  the  recognition 
differential  are  random  variables. 

Using   X(t)   to  represent  a  random  variable  corresponding  to 
an  index  time   t   and  a  subscript  to  identify  the  random 
variable,  for  a  passive  sonar  receiver,  the  signal-to-noise  ratio 
in  decibels  associated  with  a  decision  at  the  index  time  can  be 
expressed  as  follows:   XSL(t)  -  XTL(t)  -  [XNL(t)  -  XDI(t)].   In 
this  expression,   SL   represents  source  level,   TL  represents 
transmission  loss,   NL   represents  noise  level  and   DI   represen- 
ts directivity  index.   Since  signal  excess   SE   is  defined  to  be 
the  difference  in  decibels  between  the  signal-to-noise  ratio  and 
the  recognition  differential  (or  detection  threshold) ,  it  too  is 
a  random  variable  and,  for  any  decision  time   t,  one  can  write: 
(1)   XSE(t)  =  XSL(t)  -  XTL(t)  -  [XNL(t)  -  XDI(t)]  -  XRD(t). 
The  distributions  of  the  random  variables  on  the  right  side  of 
Equation  1   determine  the  distribution  of  the  signal  excess.   In 
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the  passive  sonar  detection  model  described  in  Reference  6, 
xSL(t)  '   xRd(*-)   and,  in  effect,   XNL(t)   are  normally 
distributed  random  variables  while   XTL(t)   is  a  uniformly 
distributed  random  variable.   In  the  three  signal  excess  models 
that  are  described  in  Section  VII,  all  of  the  random  variables  in 
Equation  1   are  normally  distributed. 

It  is  sometimes  convenient  to  write  Equation  1   as  follows: 

(2)  XSE(t)  =  SE(t)  +  X(t) 

In  Equation  2,   SE(t)   is  the  expected  value  of  the  signal  excess 
determined  by  the  following  expected  value  equation: 

(3)  SE(t)  =  SL(t)  -  TL(t)  -  [NL(t)  -  DI (t) ]  -  RD(t) 

where  each  term  on  the  right  represents  the  expected  value  of  the 
indicated  random  variable  and   X(t)   is  a  random  variable  that 
determines  the  stochastic  character  of  the  signal  excess.   Since 
SE(t)   is  the  mean  of   XSE(t);  by  Equation  2,  the  mean  of   X(t) 
is  equal  to  zero  and  the  standard  deviation  of   X(t)   is  equal  to 
the  standard  deviation  of   XSE(t).   If  o      represents  the 
standard  deviation  of   XSE(t)   and  the  random  variables  on  the 
right  side  of  Equation  1   are  statistically  independent,  then 

2  2  2  2  2  2 

a   =  aSL  +  aTL  +  aNL  +  aDI  +  aRD*   This  relation  has  been  used  to 
determine  a  standard  deviation  for  the  signal  excess  in 
operational  models. 
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VI.  General  Encounter  Models 

A  basic  problem  associated  with  search  modeling  is  that 
of  determining  the  probability  that  a  target  will  be  detected  by 
at  least  one  detection  system  in  an  encounter  with  one  or  more 
detection  systems  on  or  before  the  n^*1  decision  in  the  encounter 
has  occurred  regarding  the  presence  or  absence  of  a  target  in  a 
localization  region.   In  the  encounter  models  that  are  considered 
in  this  report,  during  a  search,  observations  are  made  of  a 
series  of  localization  regions.   For  localization  regions  that  do 
not  contain  a  target,  the  signal  observed  during  the  observation 
of  the  region  is  zero.   In  addition,  in  these  models,  the  time  to 
resolve  a  false  alarm  is  zero.   However,   p^   and  pf  are 
considered  to  be   determined  by  some  criterion  such  that  pf  is 
less  than  one.   Conseguently,  although  the  time  to  resolve  a 
false  alarm  is  zero,  the  cost  associated  with  a  false  alarm  is 
not  zero. 

Using  an  order  number  rather  than  a  time  to  index  a  decision 
and  a  random  variable  N  to  represent  the  decision  order  number 
at  which  detection  first  occurs,  the  cumulative  probability  of 
detection  for  an  encounter  is  given  by:   P(N  <  n)  =  S  P(N  =  i) 
where  the  sum  index   i  =  1,2,  •••  ,n.   It  also  given  by: 
P(N  <  n)  =  1  -  (l  -  gx)(l  -  g2)  •••  (1  -  gn)   where 


g^  =  P(N  =  i|N  <  i-1)   is  the  probability  of  the  event  detection 
on  the  itn  decision  conditioned  on  the  event  no  detection  on  an 
earlier  observation.   The  second  expression  for  P(N  <  n)   is 
generally  of  greater  interest  than  the  first  expression,  since 
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g^   can  usually  be  more  directly  related  to  operational 
parameters  such  as  target  range  and  target  environment  that 
determine  a  targets  detectability  than  can   P(N  =  i) .   If  a 
target  is  not  present  in  the  localization  region  that  corresponds 
to  the  itn  decision,  then  g^  will  be  zero.   The  second 
expression  above  for   P(N  <  n)   can  also  be  written  in  the 
following  form:   P(N  <  n)  =  1  -  exp[E  ln(l  -  q±) ]   where  the  sum 
index   i  =  1,2,  •••  ,m.   In  this  expression,   ln(l  -  q±)      can  be 
approximated  by   -g^   for  g^  <<  1.   With  this  approximation, 
P(N  <  n)  =  1  -  exp(-S  gi) .   This  expression  provides  a  basis  for 
developing  a  continuous  analog  for   P(N  <  n)   that  has  been  used 
to  describe  detection  systems  for  which  decisions  can  be 
considered  to  occur  continuously.   That  is,  for  cases  where  the 
time  of  the  observations  corresponding  to  a  decision  and  the  time 
between  decisions  are  both  negligible  relative  to  the  average 
time  of  an  encounter.   The  analog  can  be  developed  as  follows: 
First,  replace  the  decision  index  by  a  time  index.   Next,  let   6t 
be  the  time   between  decisions.   Then,  the  cumulative  probability 
of  detection  after  the  ntn  decision  can  be  expressed  by: 
P(N  <  n)  =  P(T  <  t)   where  T  =  N«St   and   t  =  n«6t.   Next, 
define  a  detection  rate  function  (a  probability  of  detection  per 


unit  time)  by  the  expression:   (l/6t)P[T  =  t|T  <  (n-l)»5t]. 
If   T   is  considered  to  be  a  continuous  random  variable,  then 
this  expression  can  be  used  as  a  guide  for  defining  a  continuous 
detection  rate   t (t) .   The  definition  is: 


19 


(4)  T(t)  =  lim{ (l/6t)P(t  <  T  <  t+6t|T  <  t) } 

where  the  limit  is  for   &t   approaching  zero  and   P(T  <  t)   is 
such  that   P(T  <  t)  =  P(N  <  n)   when  t   corresponds  to  n   for 
all  values  of   n  that  are  involved  and   P(T  <  t)   is  a 
nondecreasing  function  of   t.   Eguation  4   implies  the  following 
differential  equation:   dp(t)/dt  =  [1  -  p(t)]T(t)   where 

p(t)  =  P(T  <  t) .   The  solution  to  this  differential  equation  is: 

'tl 

(5)  P(T  <  tl)  =  1  -  P(T  <  tO)  exp[-Jt0  T(t)dt] 

where  t  is  the  time  index  for  a  decision  during  an  encounter,   to 
and  tl  are  times  of  the  encounter  and  tl  >  to.   A   T(t)   that  is 
based  on  a  visual  detection  model  is  described  in  Reference  7. 
If  the  detection  capability  of  a  detection  system  is  assumed  to 
depend  on  a  target's  position  relative  to  the  detection  system 
during  an  encounter  but  not  to  depend  on  the  clock  time,  then  the 
time  index  of  a  decision  can  be  a  relative  index  that  determines 
a  target  position  that  is  associated  with  a  decision  rather  than 
the  clock  time  associated  with  the  decision. 

The  above  results  apply  to  the  case  of  an  encounter  between 
a  target  and  a  detection  system  consisting  of  more  than  one 
sensor.   However,  if  the  sensors  are  not  located  on  a  single 
platform  or  at  a  single  point,  it  is  generally  more  convenient  to 
describe  encounters  of  this  kind  in  terms  of  the  encounters 
between  the  target  and  the  individual  sensors.   In  either  case, 
if  the  event  that  a  sensor  detects  a  target  during  an  encounter 
is  not  independent  of  the  this  event  for  other  sensors,  in  order 
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to  describe  this  in  an  encounter  model,  the  correlation  between 
this  sensor  and  the  other  sensors  must  be  specified.   This  has 
been  done  in  some  models  as  follows:   First  determine  the 
probability  of  detection  during  the  encounter  for  each  sensor 
involved.   Then  determine  the  probability  that  at  least  one  of 
the  sensors  will  detect  the  target  during  the  encounter  for  the 
case  where  the  random  factors  that  determine  detection  are 
independent  and  again  for  the  case  where  they  are  completely 
dependent.   Let   P-^   be  the  probability  that  the  i*-"  sensor 
detects  the  target  during  the  encounter  and   let  the  probability 
that  at  least  one  sensor  detects  the  target   be   Pj  in  the  first 
case  and   Pp  in  the  second  case.   For  the  first  case, 
Pj  =  1  -  (1  -  P^) (1  -  P2)  •••  (1  -  Pn)   where  n   is  the  number 
of  sensors  involved.   And,  for  the  second  case, 

PD  =  pm  where  Pm  >  P-^   for   i  =  1,2,  •••  ,n.   Next,  combine  the 
two  probabilities  as  follows:   P  =  a* PD  +  (1  -  a) «Pj .   A  way  to 
determine  a  value  for  a   is  described  in  Reference  8. 
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VII.   Three  Signal  Excess  Encounter  Models 

In  the  three  models  described  in  this  section, 
detection  is  defined  in  terms  of  signal  excess  as  it  is  in 
Section  V.   Each  of  the  models  determines  a  cumulative 
probability  of  detection   for  a  target  in  an  encounter  with  a 
passive  sonar  system.   An  observation  in  the  models  is  indexed  by 
time  and  the  index  can  usually  be  considered  to  be  the  time  at 
the  end  of  the  observation.   During  an  encounter,  observations 
are  made  of  a  series  of  localization  regions.   For  localization 
regions  that  do  not  contain  a  target,  the  signal  observed  during 
the  observation  of  the  region  is  zero.   The  time  to  resolve  a 
false  alarm  is  zero.   However,  false  alarms  are  not  ignored  in 
that  the  value  of  RD   (or   DT)   is  determined  by  some  specified 
false  alarm  probability.   Consequently,  although  the  time  to 
resolve  a  false  alarm  is  zero,  the  cost  associated  with  a  false 
alarm  is  not  zero.   In  determining  signal  excess  in  the  models, 
it  is  convenient  to  use  Equation  2.   For  each  decision  in  an 
encounter,  there  is  a  random  variable  X(t)   defined  by  Equation 
2   that  determines  the  random  character  of  the  signal  excess. 
For  a  sequence  of  decisions,  the  set  of  these  random  variables 
ordered  by  their  time  index  constitutes  a  stochastic  process. 
The  joint  distributions  of  these  random  variables  determines  the 
nature  of  this  process.   In  the  three  encounter  detection  models 
described  in  this  section,  the  process  is  called  a  lambda  sigma 
jump  process.    The  jumps  in  the  series  occur  at  times  determined 
by  a  Poisson  process  with  a  mean  rate  lambda.   This  implies  that 
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the  time  between  jumps  is  a  random  variable  with  an  exponential 
distribution  and  that  the  expected  times  between  jumps   t   is 
equal  to  the  reciprocal  of  lambda.   The  time  series  that  are 
generated  by  lambda  sigma  jump  processes  are  represented  by 
Figure  2   below. 


— 
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Figure  2.   A  time  series  representing  a  realization  of  a  lambda 
sigma  jump  process.   On  the  plot,   a   in   dB   equals  one  unit  on 
the  vertical  axis  and   t   equals  one  time  unit  on  the  horizontal 
axis. 

From  Figure  2,  note  that  the  observed  values  of  neighboring 
random  variables  are  equal  unless  a  jump  has  occurred  between 
them.   When  a  jump  occurs,  the  first  random  variable  after  the 
jump  is  normally  distributed  with  mean  zero  and  variance   o2   and 
it  is  independent  of  all  the  random  variables  before  the  j amp . 
Conditioned  on  a  jump  pattern,  this  random  variable  and  all  the 
random  variables  between  it  and  the  next  jump  are  dependent  and 
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the  correlation  coefficient  between  any  pair  is  one.   That  is,  if 
the  value  of  the  signal  excess  is  known  at  some  time,  then  all  of 
the  values  between  the  last  jump  before  that  time  and  the  first 
jump  after  that  time  are  also  known.  However,  since  the  jumps 
occur  randomly,  knowing  the  value  of  the  signal  excess  with 
certainty  at  some  time  does  not  determine  the  values  of  the 
signal  excess  with  certainty  at  neighboring  times.   In  the 
unconditioned  case,  the  correlation  coefficient  between  the 
random  variables  X(t)   and  X(t+T)   is  egual  to   1/e.   For  this 
reason,   t   is  referred  to  as  a  relaxation  time.   If  only  the 
transmission  loss  is  considered  to  be  determined  by  a  stochastic 
process,  then  only  the  transmission  loss  will  be  a  random 
quantity  and   t   is  the  expected  time  between  jumps  in  the  random 
component  of  the  transmission  loss.   It  appears  that  the  use  of 
the  lambda  sigma  jump  process  is  based  more  on  past  practice  than 
on  experimental  justification.   In  this  regard,  see  Reference  9. 
By  referring  to  Equation  1,  it  can  be  seen  that  the  lambda  sigma 
jump  stochastic  process  is  equal  to  the  sum  of  the  stochastic 
processes  that  determine  the  random  variables  on  the  right  side 
of  this  equation.   But,  the  nature  of  the  processes  should  depend 
on  the  encounter  that  is  being  modeled.   This  suggests  that  an 
adequate  description  of  an  encounter  using  a  lambda  sigma  jump 
process  might  not  be  possible  in  all  cases. 

In  the  three  encounter  models  described  below,  detection  is 
defined  in  terms  of  signal  excess  as  described  in  Section  IV  and 
decisions  are  indexed  by  a  time  that  can  usually  be  considered  to 
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be  the  time  of  the  decision.   During  an  encounter,  observations 
are  made  of  a  series  of  localization  regions.   For  localization 
regions  that  do  not  contain  a  target,  the  signal  observed  during 
an  observation  of  the  region  is  zero.   In  the  models,  the  time  to 
resolve  a  false  alarm  is  zero.   However,  false  alarms  are  not 
ignored  in  that  the  value  of   RD   (or  DT)   is  determined  by  some 
specified  false  alarm  probability.   Consequently,  although  the 
time  to  resolve  a  false  alarm  is  zero,  the  cost  associated  with  a 
false  alarm  is  not  zero. 

The  First  Passive  Sonar  Encounter  Detection  Model:   This 
model  describes  an  encounter  as  a  series  of  decisions.   In  the 
model,  the  integration  time  that  determines  the  recognition 
differential  is  equal  to  the  duration  of  the  observations 
corresponding  to  a  decision.    For  an  encounter  in  which  SE(t) 
is  unimodal  and  in  which  the  time  of  the  single  maximum  is  prior 
to  the  end  of  the  encounter,  it  is  shown  in  Reference  10  that 
the  probability  that  detection  will  occur  during  the  encounter 
is  given  by  the  following  equation: 

(6)    p  =  1  -  [(1  -  pc)/(l  -  p-pc)](l  -  p-Pl)  •••  (1  -  p-pm) 
where   p  =  1  -  exp(-St/T)   and  p-^  =  *[SE(ti)/o]   for   1  =  1,2,3, 
•••  ,m.   Here,   St   indicates  the  duration  of  an  observation  and 
*   indicates  the  standard  normal  cumulative  distribution  function 
as  before.   The  integer  c   is  the  index  of  the  decision  time  tc 
for  which  SE(tc)   is  greater  than  or  equal  to  SE(tjJ   for  any 
time  t^   and  tc  <  tm. 
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As   t   approaches  zero,   8   approaches  one  and  Equation  6 
approaches  this  form: 

(7)  p  =  1  -  (1  -  Pi)  ...  (1  -  pm). 

In  this  limit,  the  signal  excess  random  variables  are  all 
independent.   Note  that  Equation  7   applies  without  the  condition 
that  SE(t)   be  unimodal. 

As   t   approaches  infinity,   B   approaches  zero  and 
Equation  6   approaches  this  form: 

(8)  p  =  pc. 

In  this  limit,  the  correlation  coefficient  between  any  pair  of 
signal  excess  random  variables  is  equal  to  one.   Note  that 
Equation  8   applies  without  the  condition  that  SE(t)   be 
unimodal.   Equation  8   is  said  to  define  a  complete  dependence 
encounter  model. 

The  Second  Passive  Sonar  Encounter  Detection  Model:   This 
model  is  in  a  sense  a  third  limiting  form  of  the  first  passive 
sonar  encounter  detection  model.   In  this  limit,  the  time  between 
decisions  approaches  zero.   However,  in  this  limit  the 
integration  time  that  determines  the  recognition  differential  is 
not  equal  to   St  and  it  does  not  approach  zero.   It  is,  in 
effect,  chosen  by  the  user  of  the  model  through  the  user's  choice 
of  the  value  for  the  recognition  differential. 

For  an  encounter  that  begins  at  tl   and  ends  at  t2   and 
for  which  SE(t)   is  unimodal,  it  is  shown  in  Reference  10  that 

for  this  limit,  Equation  4   has  the  following  form: 

pt2 

(9)  p  =  1  -  [1  -  p(tc)]  exp[-(l/T)  Jtl  p(t)  dt] 
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where   p(t)  =  *[SE(t)/a]   and  where  now   tc   is  the  encounter 
time  such  that   SE(tc)   is  greater  than  equal  to   SE(t)   for  any 
other  encounter  time   t   and   t2  >  tc. 

The  Third  Passive  Sonar  Encounter  Detection  Model:    This 
model  describes  an  encounter  between  a  target  and  a  passive  sonar 
detection  system  in  which  detection  occurs  at  a  time   t   if  the 
average  value  of  the  square  of  the  continuously  observed 
signal-to-noise  ratio  over  a  time  interval  of  length  u   from 
time   t-u   to   time   t   is  greater  than  or  equal  to  the  square  of 
the  signal-to-noise  ratio  that  determines  the  recognition 
differential  for  that  time  interval.   With  R(s)   the  random 
signal-to-noise  ratio  at  a  time   s   and  Rm(u)   the  random 
signal-to-noise  ratio  that  determines  the  random  recognition 
differential  for  the  time  interval   u   (the  integration  time) , 

detection  at  time   t   is  the  following  event: 

rt 

(10)  (1/u)  Jt_u  [R(s)/R(u)m]2  ds  >  1. 

The  random  integrand  in  the  expression  that  defines  the  event  is 
related  to  the  random  signal  excess  at  the  time   s   for  an 
integration  time   u.   The  relation  is: 

(11)  10  log  [R(s)/Rm(u)]2  =  2[SE(s;u)  +  X(s)] 

where   SE(s;u)   is  the  expected  value  of  the  signal  excess  at  a 
time   s   for  an  integration  time   u   and   X(s)   is  the  random 
component  of  the  signal  excess  at  the  time   s.   In  the  model, 
X(s)   is  determined  by  a  lambda  sigma  jump  process  and   SE(s;u) 
is  determined  by  an  expected  value  sonar  equation  with  a 
recognition  differential:   RD(u)  =  10  log  rm(u) .   Here,   **m(u) 
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is  the  expected  value  of  the  signal-to-noise  ratio  that  gives  a 
probability  of  detection  equal  to   .5   for  an  integration  time  u 
and  a  specified  probability  of  false  alarm  pf.   With  the  signal 
detection  process  described  by  a  Case  II  signal  detection  model, 
the  detection  index  necessary  to  give  the  required  operating 
point   (pf,.5)   is  related  to  the  integration  time  u  and  the 
signal-to-noise  ratio   rm(u)   by: 

(12)  d  =  u-(BW) [rm(u)]* 

where   BW   is  the  bandwidth  of  the  receiver.   For  a  spectrum 
analyzer,   BW  would  be  the  bandwidth  corresponding  to  a  given 
frequency  resolution  and  d  would  be  the  detection  index 
required  in  order  to  be  at  the  operating  point   (pf,.5)   for  a 
signal  that  was  contained  within  a  bandwidth   BW.   Since  d   is 
fixed,  because  of  Equation  12: 

(13)  RD(u)  =  5  log(t0/u)  +  RD(t0) 

where   t0  is  a  reference  integration  time.   By  using  Equation  13, 
Equation  10   becomes: 

%t      (1/5) [X(s)  +  SE(s;t0)  -  5  10g(t0)] 

(14)  Jt_u  10  ds  >  1 

after  moving  the  time  origin  so  that  t  >  0.   In  this  equation, 
u  =  t   for  t  <  t0   and  u  =  t0   for  t  >  t0   and  SE(s;t0)   is 
the  expected  value  of  the  signal  excess  at  the  time   s   for  an 
integration  time  t0.   In  an  encounter,  detection  occurs  the 
first  time  that  Relation  14   is  satisfied. 

As  is  pointed  out  in  Reference  11,  the  appeal  of  the  Third 
Passive  Sonar  Encounter  Model  relative  to  the  Second  and  First 
Passive  Sonar  Encounter  Detection  Models  is  that  it  appears  to 
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more  closely  describe  the  detection  process  in  passive  sonar 
detection  systems  that  display  their  processed  data  to  an 
operator  in  a  continuous  manner.   However,  results  reported  in 
Reference  12   indicate  that  the  difference  between  the  three 
models  may  not  be  significant  in  some  types  of  encounters. 
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VIII.   Straight  Line  Encounters 

Suppose  a  target's  detectability  depends  on  its  range 
from  a  detection  system  and  that  the  probability  of  detection  is 
effectively  zero  beyond  a  range   rm   for  any  target  azimuth.   In 
this  report,  an  encounter  between  the  target  and  the  detection 
system  is  the  event  that  the  range  between  the  target  and  the 
detection  system  is  less  than  or  equal  to   rm.   In  addition, 
suppose  rm   is  small  enough  so  that  when  the  target  and  the 
detection  system  are  having  an  encounter  they  can  be  considered 
to  be  moving  on  planes  parallel  to  a  tangent  plane  to  the  earth's 
surface  at  some  point  in  their  vicinity.   If  this  is  the  case, 
then  while  the  target  and  detection  system  maintain  a  constant 
course  and  speed  during  an  encounter,  the  encounter  is  called  a 
straight  line  encounter. 

A  straight  line  encounter  can  be  described  in  terms  of  a  two 
dimensional  rectangular  coordinate  system  whose  plane  is  parallel 
to  the  tangent  plane  to  the  earth.   If  the  coordinate  system  is 
stationary  relative  to  the  detection  system  with  the  detection 
system  located  at  the  origin  and  is  oriented  so  that  the  target's 
motion  is  parallel  to  the  y-axis   and  is  in  the  positive 
y-direction,  then  the  target's   x-coordinate   during  a  straight 
line  encounter  will  be  constant.   The  constant  is  equal  to  the 
target's  horizontal  range  at  the  closest  point  of  approach  (CPA) 
on  the  straight  line  track  on  which  the  target  is  moving  relative 
to  the  detection  system  during  the  encounter.   This  range  is 
called  the  target's  lateral  range. 
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A  complete  straight  line  encounter  is  one  in  which  a  target 
begins  an  encounter  on  a  straight  track  at  a  point  on  the 
relative  track  whose  range  from  a  detection  system  is  equal  to 
rm,   continues   to   CPA  and  then  continues  to  a  point  on  the 
relative  track  whose  range  from  the  detection  system  is  again 
equal  to  rm.   Let   p(x)   be  the  cumulative  probability  that  a 
target  is  detected  by  a  detection  system  in  a  complete  straight 
line  encounter  in  which  the  target's  lateral  range  is  x.   Then 
the  function  p(x)   defines  what  is  called  a  lateral  range  curve 
or  lateral  range  function. 

Let  p  be  the  probability  that  a  target  is  detected  during  a 
complete  straight  line  encounter.   If  the  lateral  range  of  a 
target  in  a  straight  line  encounter  is  assumed  to  be  a  continuous 
random  variable   X  with  a  uniform  distribution  that  extends  at 
least  from  -rm  to   rm  where   rm   is  the  maximum  detection 
range,  then  the  probability  that  a  target  will  be  detected  during 
a  complete  straight  line  encounter  is: 

'00 

(15)      p  =  (l/2rm)  J_oo  p(x)dx 

where  the  limits  of  integration  can  be  used  since  the  value  of 
p(x)   is  zero  for   x  <  -rm  and  for   x  >  rm.      The  quantity: 


W  =   J-co  P(: 


(16)       W  =  j^o   p(x)dx 

is  called  the  sweep  width.   Using  the  definition  of  W  given  by 
Equation  16,  the  probability  that  a  target  will  be  detected 
during  a  complete  straight  line  encounter  under  the  conditions 
that  determine  Equation  15   is   p  =  (l/2rra) *W. 
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IX.  Two  Intermittent  Signal  Encounter  Models 

In  the  intermittent  signal  encounter  models  that  are 
described  in  this  section,  an  encounter  is  a  complete  straight 
line  encounter  and  during  an  encounter  a  target  emits  a  signal  at 
various  times.   Two  cases  are  considered:   In  the  first  case,  the 
signals  are  emitted  periodically,  the  signals  are  of  length   St 
and  the  time  between  signals  is   t   where   t  >  St.   In  the  second 
case,  the  signals  are  instantaneous  and  the  signals  are  emitted 
at  times  determined  by  a  Poisson  process  for  which  the  expected 
time  between  signals  is  equal  to   t.   In  the  model,  the 
detectability  of  a  target  signal  depends  on  a  target's  horizontal 
range  from  a  detection  system,  but  on  no  other  factors.   If  a 
signal  is  emitted  while  the  target  is  within  a  range   r,  it  will 
be  detected.   For  a  continuous  signal,  the  lateral  range  function 
of  a  detection  system  for  a  target  is:   p(x)  =  1   for   |x|  <  r 
and   p(x)  =  0   for   |x|  >  r  where  the  horizontal  range   r   is 
determined  by  the  characteristics  of  the  detection  system  and  the 
target.   The  geometry  for  an  encounter  is  shown  in  Figure  1 
below. 

For  an  intermittent  signal,  the  length  of  a  target's  track 
relative  to  a  detection  system  where  if  a  signal  is  emitted  it 
will  be  detected  is   2(r*  -  x1)*  +  wSt  with  w  the  speed  of 
the  target  relative  to  the  detection  system.   So,  a  target's 
exposure  time  during  an  encounter  is   (2/w) (r*  -  x2)*  +  St. 

For  a  periodic  signal,  there  are  two  cases.   In  the  first 
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case,   r  >  w(t  -  St)/2.   In  this  case,  the  exposure  results  in 
the  following  lateral  range  function: 

p(x)  =  0   for   |x|  >  r 
(17)       p(x)  =  1   for   |x|  <  (rJ  -  [w(t  -  St)/2]1}^ 

p(x)  =  (2/wT)(r2  -  x*)2  +  5t/T   otherwise 

In  the  second  case,   r  <   w(t  -  6t)/2.   In  this  case,  the  middle 

eguality  in  Eguation  17   does  not  apply. 

For  a  signal  that  is  instantaneous  and  whose  emission  times 

are  determined  by  a  Poisson  process,  the  exposure  results  in  the 

following  lateral  range  function: 

p(x)  =  1  -  exp{-[ (2/w-t) (r*  -  x1)*]}   for   |x|  <  r 
(18) 

p(x)  =  0   for   |x|  >  r 

for  all  values  of   r.   A  case  in  which  the  signal  length   &t   is 

greater  than  zero  is  discussed  below. 


Figure  3.   The  encounter  geometry  for  the  two  intermittent  signal 
models  described  above. 
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For  signals  whose  emission  times  are  determined  by  a  Poisson 
process  and  whose  length  is  &t,  signals  can  overlap.   If  this  is 
allowed,  then  Eguation  18   can  be  modified  to  describe  this  case 
by  adding   5t/r   to  the  term  in  the  exponent  of  Equation  18   that 
is  within  the  square  brackets.   In  particular,  note  that  this 
modified  Equation  18   can  be  approximated  by  the  bottom  equality 
in  Equation  17   when   (2/w«T)(r*  -  x*)5  +  5t/T  <<  1.   This 
implies  that  when  the  expected  time   t   between  signals  is  large 
relative  to  the  exposure  time   (2/w) (r*  -  x2)5  +  6t  the  periodic 
emission  model  and  the  random  emission  model  are  effectively 
equivalent. 
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X.   A  Random  Search  Model 

A  search  of  a  region  in  which  a  target  moves  on  a  track 
that  consists  of  a  number  of  straight  segments  placed  in  such  a 
way  that  in  a  limiting  sense  every  section  of  the  region  is 
equally  likely  to  be  searched  on  a  segment  is  referred  to  as  a 
random  search  in  Reference  6.   Representation  of  a  search  region 
with  the  track  segments  that  could  be  imagined  to  be  the  tracks 
of  a  random  search  are  shown  in  Figure  4 . 


Figure  4.   A  search  region  and  a  track  that  could  be  described  as 
a  random  search  track. 

Two  developments  of  a  model  to  describe  this  kind  of  search 
are  contained  in  this  section.   The  first  development  is  based  on 
the  following  conditions:   1.  A  target  is  at  a  fixed  position 
within  a  defined  search  region.   2.  A  searcher's  track  is  a 
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sequence  of  straight  line  segments  that  are  within  the  search 
region.   3.  The  searcher's  detection  system  is  such  that  while  on 
a  track  segment,  a  rectangle  is  searched  that  is  contained  within 
the  search  region,  is  of  length  equal  to  the  length  of  the  track 
segment  and  is  oriented  so  that  its  long  axis  is  parallel  to  the 
track  segment.   4.  The  probability  that  the  searcher's  detection 
system  will  detect  a  target  while  on  a  track  segment  with  a 
search  rectangle  that  does  not  include  the  target  is  zero.   The 
probability  that  the  searcher's  detection  system  will  detect  a 
target  while  on  a  track  segment  with  a  search  rectangle  that 
includes  the  target  is  p(x)   where  x  is  the  target's  lateral 
range  for  the  track  segment  and  p(x)   is  the  lateral  range  curve 
for  a  complete  straight  line  encounter  with  this  lateral  range.  A 
representation  of  a  search  rectangle  is  shown  in  Figure  5  below. 
5.  The  track  segments  are  located  in  such  a  way  that  the  event 
that  the  target  is  within  the  search  rectangle  associated  with  a 
track  segment  is  independent  of  the  event  that  the  target  in  the 
search  rectangle  associated  with  any  other  track  segment.   And 
the  probability  of  the  event  is  equal  to  the  ratio  of  the  area  of 
the  search  rectangle  to  the  area  of  the  search  region  and,  given 
a  target  is  within  a  search  rectangle,  its  position  is  uniformly 
distributed  over  the  rectangle. 

Condition  4   implies  that  the  random  search  model  is  based 
on  the  concept  of  a  complete  straight  line  encounter.   The 
definition  of  an  encounter  that  is  intended  here  is  that  given  in 
Section  VI.   This  implies  that  in  the  random  search  model  the 
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time  to  resolve  a  false  alarm  is  zero.   However,  for  the  model, 
P(j   and  pf  are  considered  to  be  determined  by  some  criterion 
such  that   pf  is  less  than  one.   Consequently,  although  the  time 
to  resolve  a  false  alarm  is  zero  in  the  model,  the  cost 
associated  with  a  false  alarm  is  not  zero.   (A  concise  model  that 
accounts  for  the  time  to  resolve  false  alarms  is  described  in 
Reference  13.)   Condition  4   also  implies  that  when  a  searcher  is 
on  a  track  segment  with  a  search  rectangle  that  contains  a 
target,  the  encounter  is  a  complete  straight  line  encounter.   And 
Condition  5,  which  can  be  considered  to  specify  a  random 
arrangement  of  the  track  segments,  implies  that  when  this  is  the 
case,  for  the  complete  straight  line  encounter,  the  target's 
lateral  range  is  a  random  variable  that  is  uniformly  distributed 
between   -b/2   and  b/2   where   b   is  width  of  the  search 
rectangle  (the  dimension  of  the  rectangle  perpendicular  to  the 
associated  track  segment) . 


Swept  Area 


Track 


Figure  5.   A  track  segment  and  its  associated  search  rectangle 
that  could  correspond  to  a  search  with  an  aircraft  mounted 
infrared  detection  system. 
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Based  on  the  above  considerations,  the  probability  that  a 
target  will  be  detected  while  a  searcher  is  on  a  track  segment 
with  an  associated  search  rectangle  that  contains  the  target  is 
given  by: 

(19)  J_oo  p(x)fx(x)dx  =  W/b 

where   fx(*)  =  1/2   for   -b/2  <  x  <  b/2   and   fx(*)  =  °   and 
p(x)  =  0   otherwise.   Note  that  the  left  side  of  Equation  19 
applies  to  any  complete  straight  line  encounter  in  which  the 
target's  lateral  range  for  the  encounter  is  considered  to  be  a 
random  variable  with  a  distribution  determined  by  the  probability 
density  function   fx(x)  •   If  it  is  not  given  that  the  target  is 
within  the  search  rectangle  associated  with  a  track  segment,  then 
the  unconditional  probability  that  the  target  will  be  detected  on 
the  track  segment  is  given  by:   (W/b) (SA/A)   where   &A   is  the 
area  of  the  search  rectangle  associated  with  the  track  segment 
and  A  is  the  area  of  the  search  region.   With   1   the  length  of 
the  rectangle,   SA  =  b»l   and  the  probability  becomes:   (W»l)/A. 
Then,  since  the  event  that  the  target  will  be  in  the  search 
rectangle  of  a  track  segment  is  independent  of  the  event  that  it 
will  be  in  the  search  rectangle  of  any  other  track  segment,  the 
probability   p   that  a  random  search  consisting  of  m  track 
segments  will  detect  the  target  is  given  by: 

1  -  [1  -  (W«l1)/A][l  -  (W'l2)/A]  •••  [1  -  (W-ln)/A]   where   li 
is  the  length  of  the   itn  track  segment.   The  probability  is 
also  given  by:   p  =  1  -  exp{S  ln[l  -  (W*  lj_)/A]  }  .   If  for 
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i  =  1,2,  •••  ,n   (W'l^)/A  <<  1,   then  this  expression  can  be 

approximated  by: 

(20)       p  =   1  -  exp[-(W»l)/A] 

where  the  sum  index   i  =  1,2,  •••  ,m  and   1=21^   is  the  track 

length  of  the  search.   Equation  2  0   is  known  as  the  random  search 

formula. 

The  second  development  is  based  on  the  following  condition: 
The  detection  rate  of  a  searcher's  detection  system  for  a  target 
is  given  by:   t (t)  =  [W«v(t)]/A  where   t (t)   is  the  continuous 
detection  rate  function  in  Equation  5.   Using  the  value  for   t (t) 
and  Equation  5,  the  random  search  formula  given  by  Equation  18  is 
easily  found. 

The  advantage  of  the  method  used  in  the  second  development 
of  the  random  search  formula  over  the  method  used  in  the  first 
development  is  that,  in  general,  it  can  be  more  easily  used  to 
determine  the  probability  of  detection  for  a  search  that  can  be 
defined  in  terms  of  a  detection  rate.   Reference  14   contains  an 
example  of  an  application  of  Equation  5   to  a  case  where  the 
search  region  expands  with  time. 


39 


XI.   Ladder  and  Barrier  Search  Models 

In  some  barrier  searches,  the  barrier  search  track  is  a 
ladder  search  track  relative  to  a  reference  system  that  moves 
with  the  target.   This  fact  is  used  in  the  barrier  search  model 
development  that  follows  the  two  ladder  search  model  developments 
below.   The  first  ladder  search  model  is  referred  to  as  an  ideal 
ladder  search  model  because  of  the  idealizations  that  are 
involved  in  its  description  of  a  ladder  search.   The  second 
ladder  search  model  is  referred  to  as  a  degraded  ladder  search. 
It  can  be  considered  to  describe  a  ladder  search  track  in  which 
navigational  errors  result  in  omissions  and  overlaps  in  coverage. 

An  Ideal  Ladder  Search  Model:   The  model  is  based  on  the 
following  conditions:   1.  A  ladder  search  region  is  a  rectangle 
that  contains  a  fixed  target.   2.  During  a  search  of  the  region, 
the  searcher's  detection  system  moves  on  a  set  of  m  parallel 
track  segments  of  length  b   separated  by  a  distance   s.   3 .  As 
the  detection  system  moves  along  a  track  segment,  it  searches  a 
rectangular  strip  of  length  b   and  width   s   within  the  search 
region.   4.  The  m   rectangular  strips  that  correspond  to  the  m 
track  segments  completely  cover  the  ladder  search  region  with  no 
overlap.   5.  If  a  target  is  within  the  rectangular  strip 
corresponding  to  a  track  segment,  then  there  will  be  a  complete 
straight  line  encounter  between  the  target  and  the  detection 
system  when  the  detection  system  moves  along  the  track  segment 
and  the  lateral  range  of  the  encounter  will  be  uniformly 
distributed  across  the  width  of  the  strip.   If  the  target  is  not 
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in  the  rectangular  strip,  then  there  will  not  be  an  encounter  and 
the  probability  that  the  target  will  be  detected  while  the 
detection  system  is  on  the  track  segment  is  zero. 


—  s  — 

Figure  6.   A  schematic  representation  of  a  ladder  search  geometry 
for  a  case  in  which  the  ladder  search  track  segments  are 
superimposed  on  and  bisect  their  corresponding  rectangular 
strips. 


Since  targets  outside  of  the  rectangular  strip  that 
corresponds  to  a  track  segment  cannot  be  detected  while  a 
detection  system  is  on  the  track  segment  because  of  Condition  5, 
in  the  model,  the  sweep  width  W   of  a  searcher's  detection 
system  must  satisfy  the  relation   W  <  s  .   In  particular,   W  =  s 
only  holds  when  the  detection  system  detects  a  target  that  is  in 
a  rectangular  strip  with  probability  one  for  any  target  lateral 
range.   This  kind  of  detection  system  is  sometimes  referred  to  as 
a  cookie  cutter  detection  system.    However,  this  terminology  can 
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be  misleading  since  it  suggests  the  detection  system  detects 
equally  well  for  all  azimuths.   But  this  is  not  a  requirement  on 
the  system  in  order  that  W  =  s. 

The  ideal  ladder  search  model  implies  that  if  the  conditions 
of  the  model  are  satisfied,  then  the  probability  p   that  a 
target  will  be  detected  by  a  an  ideal  ladder  search  is  given  by: 

(21)  p  =  W/s 

where  W/s   <  1.   The  quantity  W/s   is  called  the  coverage 
factor  in  this  case. 

A  Degraded  Ladder  Search  Model:      The  above  model  implies 
perfect  navigation  in  addition  to  other  idealizations.   A  model 
of  a  ladder  search  is  given  in  Reference  6   that  could  be  used 
for  cases  in  which  this  is  a  poor  assumption.   The  model  which  is 
referred  to  here  as  a  degraded  ladder  search  model  can  be 
considered  to  describe  navigational  inaccuracies  in  terms  of 
omissions  and  overlaps  of  the  rectangular  strips.   It  can  be 
developed  as  follows:   Consider  a  random  search  in  the  ladder 
search  region  whose  track  length  is  equal  to  the  search  track 
length  required  to  complete  an  ideal  ladder  search,  that  is,  a 
track  length   1  =  m*b.   The  degraded  ladder  search  model 
describes  the  result  of  omissions  and  overlaps  in  a  ladder  search 
to  be  such  that  the  probability  of  detection  for  this  random 
search  is  equal  to  the  probability  of  detection  for  the  degraded 
ladder  search.   Consequently,  since  the  area  of  the  ladder  search 
region  is  m»s«b,   for  the  degraded  ladder  search  model: 

(22)  p  =  1  -  exp(-W/s) . 
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Here,  the  requirement  that  the  coverage  factor   W/s  <  1   for 
Equation  21   can  be  relaxed.   However,  it  should  still  be 
considered  as  an  approximate  condition. 

The  condition  that  the  target  be  fixed  within  the 
rectangular  search  region  is  critical  to  both  Equation  21   and 
Equation  22.   However,  these  results  are  also  applicable  to  a 
search  for  a  moving  target  under  the  conditions  that  are 
described  next. 

A  Barrier  Search  Model:   A  target  moves  with  a  constant 
course  and  a  constant  speed   u.   Both  the  target's  course  and  the 
target's  speed  are  known  by  a  searcher.   The  searcher  establishes 
a  barrier  of  width  b   that  is  perpendicular  to  the  target's 
track  and  moves  on  the  barrier  with  a  speed  v  >  u.   The  barrier 
is  designed  so  that  in  a  reference  system  relative  to  the  target 
the  barrier  search  is  a  ladder  search  that  satisfies  the 
conditions  for  a  ladder  search  that  are  given  above.   There  are 
two  cases  to  consider:   1.  The  barrier  is  established  in  front  of 
the  target.   2.  The  barrier  is  established  behind  the  target. 

From  the  search  geometry  for  a  barrier  established  in  front 
of  the  target,  it  can  be  seen  from  Figure  7  below  that 
9  =  sin_1(u/v)  and  d  =  vt  where  t  =  s/(v  +  u)  is  the  time  to 
move  from  one  search  leg  to  the  next.  The  angle  0  and  the 
perpendicular  distance  d  that  depend  on  u,  v  and  s,  and  the 
width  of  the  barrier  b  are  the  quantities  that  are  required  in 
order  to  establish  the  barrier  operationally. 
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Searcher 


Figure  7.  A  barrier  search  track  shown  for  a  barrier  established 
in  front  of  the  target.  The  track  is  shown  in  a  reference  system 
fixed  relative  to  the  earth. 

For  a  barrier  that  is  established  in  front  of  a  target,  one 
of  three  barrier  types  will  result.   A  barrier's  type  is 
determined  by  the  relation  of  the  distance  d  to  the  distance 
g  =  ut  where  the  time  t  =  b/(v*  -  u*)5   is  the  time  to  complete 
a  search  leg  (cross  the  barrier) .   The  barrier  type  is  determined 
as  follows:   1.  For  g  <  d,  the  barrier  is  an  advancing  barrier. 

2.  For  g  =  d,  the  barrier  is  a  stationary  barrier. 

3.  For  g  >  d,  the  barrier  is  a  retreating  barrier. 

For  a  barrier  established  behind  the  target,  there  is  only 
one  barrier  type  and  it  is  called  an  overtaking  barrier.   For  an 
overtaking  barrier,   6  =  sin~1(u/v)   as  for  a  barrier  established 
in  front  of  the  target.   But,  for  an  overtaking  barrier, 
t  =  s/(v  -  u)   and  d  =  vs/(v  -  u)  . 


44 


Given  the  target  crosses  the  barrier,  the  probability  of 
detection  for  an  ideal  barrier  search  is  given  by  Equation  21 
and  the  probability  for  a  degraded  barrier  search  is  given  by 
Equation  22  where  the  terminology  refers  to  the  nature  of  the 
ladder  search  in  the  reference  system  moving  with  the  target. 
A  discussion  of  an  application  of  these  two  equations  to  a  search 
for  a  magnetic  anomaly  target  is  given  in  Reference  15. 
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XII.   A  Target  State  Estimation  Procedure 

A  target  state  estimation  procedure  based  on  bearing 
observations  is  developed  in  this  section  that  generates  point 
estimates  of  a  target's  position  and  velocity  vector  coordinates 
in  a  rectangular  coordinate  system.   The  procedure  is  based  on  a 
model  in  which  bearing  errors  are  unknown  and  are  not  determined 
by  random  variables  with  known  distributions.   Because  of  this, 
confidence  regions  for  the  estimates  are  not  generated  by  the 
procedure.   However,  for  a  moving  target,  it  illustrates  some 
general  characteristics  of  bearings  only  target  motion  analysis 
(TMA) . 

The  target  state  estimation  procedure  model  is  defined  by: 
1.  The  target  moves  in  a  plane  with  a  constant  but  unknown  course 
and  speed.   2.  Observations  of  the  target  are  made  from  known 
positions  at  known  times.   3.  The  observations  provide  only 
target  bearings  with  unknown  errors.   4.  Target  position 
estimates   [xt (i) ,yt (i) ]   for   i  =  1,2,  •••  ,  n   and  target 
velocity  component  estimates  ux  and  uy  are  chosen  so  that  the 
sum  S  =  S  d{      is  a  minimum  where  d^   is  the  algebraic  distance 
between  the  estimated  position  and  the  bearing  line  corresponding 
to  the  observed  bearing. 

Because  of  the  reguirement  that  the  target  move  with 
constant  course  and  speed  during  the  encounter,  the  number  of 
independent  estimates  is  reduced  from   2m  to   4:   Xfc(l),  Yt^1) ' 
ux  and  Uy.   These  four  estimates  constitute  a  set  of  target 
state  parameters  that  define  a  target's  motion  by  the  relations: 
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xt(i)  =  xt(x)  +  ux'(ti  "  tl)      and  yt(i)  =  Yt(1)  +  Uy-(ti  -  tx) . 
The  geometry  on  which  the  procedure  is  based  is  shown  in  Figure  8 
where  di  =  [xt(i)  -  x0(i)]-cos  Qi   -  [yt(i)  -  y0(i)]-sin  9^. 

North 


observed   6 
bearing 


d^  =  r^-sin  (9^  -  9^) 

[xt(i),yt(i)] 
r^   estimated  range 
9^   estimated  bearing 


[x0(i) fy0(i) ]    observer 


East 


Figure  8.   The  geometry  of  the  target  motion  analysis  model 


To  determine  estimates  of  the  target  state  parameters,  take 
the  partial  derivative  of   S   with  respect  to  each  of  them.   Then 
set  the  four  partial  derivatives  equal  to  zero.   This  creates 
four  linear  equations  in  the  four  unknown  estimates.   In  matrix 
notation,  the  equations  are:   AX  =  B  where  the  elements  of   X 
are   x-t-(l)  ,   yt(l)»   ux>   and  uy   A  necessary  condition  for  a 
unique  solution  for   X   is  that  m  >  4.   Otherwise,  the 
determinant  of  A  will  be  equal  to  zero. 

Now,  suppose  the  observations  are  at  positions  and  times 
that  correspond  to  the  positions  and  times  of  an  observer  moving 
on  some  constant  course  at  some  constant  speed  (including  zero 
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speed) .   In  this  case,  the  observation  position  coordinates  are 
related  by  the  following  equations:   x0(i)  =  x0(l)  +  vx(t-^  -  t]_) 
and  y0(i)  =  y0(l)    +  vy(t^  -  tx)   where  vx   and  vy  are  the 
required  velocity  components  of  the  observer.   Using  these 
equations  of  motion,  the  matrix  equation  AX  =  B  can  be 
transformed  to  the  matrix  equation  AX  =  0   where   the  elements 
x^2   of  the  matrix  X  are  related  to  the  elements  of  the  matrix 
a'  by  the  equations:   x^  =  xt(l)  -  x0(l)  ,   X21  =  ytU)  "  Yo(1)' 

x31  =  ux  ""  vx  and  x41  =  uy  ~  vy 

Since  the  equations  represented  by  AX  =0   are  homogenous, 

the  equations  do  not  have  unique  solutions  and  consequently  the 

equations  represented  by  AX  =  B  do  not  either.   However,  if 

there  is  at  least  one  observation  whose  time  and  position  is  not 

determined  by  the  above  equations  of  motion,  then  the 

transformation   from  X   to   X   cannot  be  made,  and  in  general  a 

unique  solution   for  X   can  be  found.   If  the  observations  are 

made  from  a  platform  that  is  moving  with  a  constant  course  and 

speed,  this  condition  can  be  achieved  by  either  changing  the 

course,  the  speed  or  both  prior  to  completing  the  observations. 

Estimation  models  that  describe  bearing  error  in  terms  of  a 

random  variable  with  a  specified  distribution  provide  a  basis  for 

determining  a  confidence  region  to  associate  with  the  point 

estimates.   A  model  is  developed  in  Reference  16  that  does  this, 

however,  it  requires  that  either  observations  of  a  target  bearing 

be  made  from  two  or  more  locations  simultaneously  or  that  the 

target  be  stationary  relative  to  the  observation  points. 

48 


XIII.   Position  Distributions  That  Change  with  Motion 

In  this  section,  target  motion  models  are  discussed 
that  provide  a  basis  for  determining  position  distributions  that 
change  with  time.   In  the  first  model,  a  target  moves  in  a  plane 
with  a  constant  course  and  speed.   In  the  second  model,  a  target 
moves  in  a  plane  with  a  constant  course  and  speed  during  a  time 
step  and  at  the  end  of  the  time  step  changes  either  its  course  or 
speed  or  both. 

Position  Distributions  based  on  the  First  Motion  Model:   Two 
cases  of  this  model  are  considered  that  provide  the  basis  for 
determining  a  position  distribution  analytically.   In  both  cases, 
the  target's  position  coordinates  in  a  rectangular  coordinate 
system  are  determined  by  a  circular  normal  distribution  at  time 
zero.   And,  in  both  cases,  during  the  motion  the  target's 
velocity  does  not  change. 

The  First  Case:   In  this  case,  the  targets 's  coordinates 
X(t)   and   Y(t)   at  time   t   are  determined  by  the  vector  sum 
with  components:   X(t)  =  X(0)  +  Ux«t   and   Y(t)  =  Y(0)  +  Uy«t 
where   X(0)   and   Y(0)   are  the  target's  coordinates  at  time  zero 
and  Ux   and  Uy   are  the  target's  velocity  components.   The 
coordinates   X(0)   and   Y(0)   are  determined  by  a  circular  normal 
distribution  whose  mean  vector  is   (0,0)   and  whose  standard 
deviation  is  a.      And  the  velocity   components   Ux   and  Uy   are 
determined  by  a  circular  normal  distribution  whose  mean  vector  is 
(ux,Uy)   and  whose  standard  deviation  is   au.   The  random  vectors 
[X(0),Y(0)]   and   (Ux»t,Uy»t)   are  independent  and  both  are 

49 


normal  random  vectors.   Since  they  are  normal  random  vectors,  the 
random  vector   [X(t),Y(t)]   is  also  a  normal  random  vector.   And, 
since  they  are  independent,  its  characteristic  function  is  equal 
to  the  product  of  their  characteristic  functions.   This  implies 
that   [X(t),Y(t)]   has  a  circular  normal  distribution  whose  mean 
vector  is   (ux«t,uv't)   and  whose  variance  is   (a*  +  a^) .   For  a 
more  complete  discussion  of  the  basis  for  the  above  argument,  see 
Reference  7. 

The  Second  Case:   The  procedure  that  is  used  to  determine 
the  distribution  of  the  random  vector   [X(t),Y(t)]   for  the 
second  case  can  also  be  used  to  determine  it  for  the  first  case. 
It  is  based  on  the  following  consideration:   In  general,  the 
joint  density  function  of  a  target's  coordinates  at  some  time  t 
greater  than  a  zero  of  time  is  determined  as  follows: 


00  poo 


(23)   fX(t)  ,  (t)  (*,y;t)  =  J-ooJ-oofx(O)  ,Y(0)  (r,s;0) fUX/Uy(v,w)dvdw 
where   r  =  x  -  vt   and   s  =  y  -  w*t  with  v  =  ux  and  w  =  uy. 
Equation  23   can  be  interpreted  as  follows:   After  multiplying 
the  right-hand  side  of  Equation  23   by   SrSs,  in  the  limit,  it  is 
the  sum  over  all  pairs  of  values   (v,w)   of  the  product  of  the 
probability  that  the  target's  velocity  components   (v,w)   are  in 
an  element   6v5w  and  the  probability  that  the  target's  position 
coordinates   (r,s)   at  time  zero  are  in  an  element  of  area   SrSs. 
This  is  the  position  such  that  the  targets' s  position  coordinates 
(x,y)   at  time  t  will  be  in  an  element  of  area   6x5y  given 
that  the  target's  velocity  is   (v,w) .   And,  in  the  limit,  this  is 
equal  to  the  left  side  of  Equation  23   after  multiplying  the  left 
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side  by   SxSy.   But,  since   5r&s   is  equal  to   SxSy   in  the 
limit,   Equation  23   results.   In  effect,   SrSs   is  translated 
from   (r,s)   to   (x,y)   without  rotation  or  distortion. 

In  the  first  case,  the  joint  distribution  of  the  of  the 
velocity  components   (Ux,Uy)   is  circular  normal  with  mean  vector 
(ux,Uy)   and  standard  deviation   au  and  the  joint  distribution 
of   [X(0),Y(0)]   is  circular  normal  with  mean  vector   (0,0)   and 
standard  deviation   o.   In  the  second  case,  the  target's  speed 
u   is  known  but  its  course  is  equally  likely  to  have  any  value 
between   0  and   2ir,  so  Ux  =  u*sin  *   and  Uy  =  u»cos  *   where   * 
is  a  uniform  random  variable  and   0  <  *  <  2tt.   In  this  case,  only 
the  single  random  variable   *   is  required  to  determine  the 
target's  velocity.   The  integral  of  Equation  23   in  polar 
coordinates  is  a  single  integral  over   #   and  the  integrand  is 
(l/2iro*)  'exp[-(r*  +  s*)/2o*  ]  •  (l/2ir)    where   r  =  x  -  u«t*sin  <f> 
and   s  =  y  -  u»t«cos  <f>.      The  result  of  the  integration  is: 
(24)   (1/2W01) •exp(-[x2  +  y*  +  (u*t) 2]/2o*)Io[ (x*  +  y2)*U't/c*] 
where   I0   indicates  the  hyperbolic  Bessel  function  of  zeroth 
order.   Expression  24   is  plotted  in  Reference  6   for  several 
values  of   t.   There,  a  single  variable  is  substituted  for  the 
distance  of  the  target  from  the  origin   (x*  +  y*)'  .   It  is 
measured  in  units  of   a   and  the  time  is  measured  in  units  of 
o/u.   The  plot  of  the  expression  shows  a  limiting  characteristic 
of  the  distribution  that  can  be  indicated  as  follows:   Multiply 
and  divide  the  resulting  expression  by  exp(-z 'U»t/a* )   where 
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z  =  (x*  +  y*)*,  the  target's  range  from  the  origin.   The  result 
of  the  operation  can  be  expressed  in  the  following  form: 

(25)  l/(2ira1)exp{-[l/(2a1)  ]  (z  -  u«t)  •  }I0  (z  'U't/a2 )  exp(-r«u«t/a2 ) 
where  the  second  factor  in  Expression  25   is  proportional  to  the 
density  function  of  a  normal  random  variable   Z   whose  mean  is 
u»t   and  whose  standard  deviation  is   a.   The  last  two  factors 
are  such  that  their  product  changes  slowly  relative  to  the  second 
factor.   The  consequence  of  this  is  that  a  plot  of  Expression  25 
against   z   in  units  of   a   for  values  of   t   in  units  of   a/u 
has  the  appearance  of  the  plot  of  a  normal  density  function  for 
values  of   t  greater  than  four  in  units  of   a/u.   In  particular, 
this  indicates  that   u*t   in  the  exponent  of  the  second  factor  in 
Expression  25   can  be  viewed  as  the  radius  of  an  average  furthest 
on  circle. 

A  target's  rectangular  coordinates  and  its  range  and  bearing 
from  the  origin  are  related  by:   X  =  R  sin  9   and   Y  =  R  cos  9. 
By  using  these  equations,  Expression  24   can  be  transformed  to: 

(26)  (1/2tt)  (r/2Trol)exp{-[r*  +  (ut)  *  ]/2a2  )  I0  (rut/a* ) 

which  is  the  joint  density  function  of  the  random  variables  R 
and  9.   The  marginal  density  function  of   R  can  be  obtained 
from  Expression  26   by  integrating  this  joint  density  function 
over  the  possible  values  of  9  which  in  radians  is  over  the 
interval   0   to  2tt.   Consequently,  the  marginal  density  for  R 
can  be  obtained  from  Expression  26  by  multiplying  it  by  2ir. 
Tabulated  values  of  the  cumulative  distribution  function  FR(r;t) 


52 


for  the  marginal  distribution  of   R   for  the  second  case  are 
listed  in  Reference  17. 

Position  Distributions  Based  on  the  Second  Motion  Model:   In 

the  second  motion  model  a  target's  motion  is  not  restricted  to  a 
constant  course  and  speed  throughout  the  motion.   At  the  end  of  a 
time  step,  a  target's  course  or  speed  or  both  can  change.   In 
general,  this  leads  to  a  position  distribution  that  can  not  be 
described  analytically.   And,  in  general,  a  monte  carlo 
simulation  method  is  required  in  order  to  determine  a  position 
distribution  that  is  based  on  such  a  model.   As  an  example  of  a 
case  in  which  the  distribution  can  be  described  analytically,  see 
Reference  18.   As  an  example  of  a  case  in  which  it  can  not, 
suppose  a  target's  initial  position  is  described  in  terms  of  a 
number  assigned  to  a  subregion  in  the  xy-plane  where  the  number 
assigned  represents  the  probability  that  the  subregion  contains 
the  target  at  an  initial  time.   In  addition,  suppose  for  each 
subregion  a  course  and  speed  distribution  is  determined  by 
assigning  numbers  to  course  and  speed  pairs  where  a  number 
represents  the  probability  the  target  will  have  the  course  and 
speed  at  the  initial  time  given  it  is  in  the  subregion  at  that 
time.   Next  suppose  for  each  course  an  speed  pair  there  is  a  time 
distribution  that  determines  the  duration  of  the  course  and  speed 
pair  and  that  the  time  distribution  is  determined  by  a  number 
assigned  to  each  discrete  time  point  where  the  number  represents 
the  probability  that  the  targets  course  and  speed  pair  will  be 
determined  by  a  new  course  and  speed  distribution  given  it  has 
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not  been  determined  at  an  earlier  time  point.   By  extending  this 
kind  of  procedure  and  then  implementing  it  in  a  monte  carlo 
simulation,  one  can  generate  complex  position  distributions  that 
describe  a  target's  position  at  discrete  time  points  in  terms  of 
probabilities  assigned  to  subregions  of  a  region  in  which  the 
target. 
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XIV.   Position  Distributions  That  Change  with  Search 

Suppose  a  target's  position  at  some  time  is  described 
by  a  position  distribution.   Now  suppose  information  becomes 
available  that  a  search  has  been  conducted  for  the  target  and 
that  the  target  has  not  been  detected.   Or  suppose  the 
information  is  that  the  target  has  been  detected.   In  the  first 
case,  negative  information  is  available  that  can  be  used  to 
modify  the  position  distribution.   In  the  second  case,  positive 
information  is  available  that  can  be  used  to  modify  the  position 
distribution.   In  both  of  the  cases,  a  target's  position 
distribution  is  assumed  to  be  specified  in  terms  of  a  set  of 
discrete  probabilities  where  each  probability  corresponds  to  a 
subregion  of  the  region  that  contains  the  target  and  each  is  the 
probability  that  the  target  is  in  that  subregion. 

Position  Distributions  and  Negative  Information:   For  a 
region  that  contains  a  target  and  consists  of   n   subregions,  let 
the   event   S^  =  {the  target  is  in  the  itn  subregion}.   And  let 
the  event   C  =  {no  detection  during  a  search  of  the  region). 
Then,  given  an  unsuccessful  search  in  the  region,  the  targets 's 
position  distribution  can  be  modified  as  follows: 
(26)    P(Si|C)   =   P(C|Si)P(Si)/P(C) 

where   i  =  1,2,  •••  ,n   and   P(C)  =  Z  P(C| Sj ) P(Sj )   with  the  sum 
from   j  =  1,2,  •••  ,n.   Note  that  Equation  26   can  be  obtained  by 
using  Bayes  theorem.   To  illustrate  how  Equation  26  might  be 
used,  suppose  that  a  search  in  a  subregion  is  considered  to  be  a 
random  search  and  that  the  sweep  width  of  detection  system 
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against  the  target  depends  on  subregion  being  searched.   For  this 
case,  let  A^   be  the  area  of  the  itn  subregion  and  let  W^   be 
the  sweep  width  in  that  subregion.   Then,  given  an  unsuccessful 
search,   P(c|Si)  =  exp(-Wj/  lj/A^)   where  1±      is  the  track  length 
of  the  searcher  in  the   ith   subregion.   Given  values  for  P(S^) , 
&i   and  1^   for   i  =  1,2,  •••  ,n   a  position  distribution  can  be 
determined  that  has  been  modified  by  the  negative  information. 

Position  Distributions  and  Positive  Information:   As  above, 
let  S^   the  event  that  a  target  is  in  the  it"  subregion  of  a 
region  that  contains  the  target.   However,  in  the  case  of 
positive  information,  instead  of  the  event  C  occurring,  the 
event   C  =  {a  contact}   occurs.   This  event  is  the  union  of  two 
mutually  exclusive  and  exhaustive  events:   T  =  {a  true  contact} 
and  F  =  {a  false  contact}.   Relative  to  the  Venn  diagram  of 
Figure  1,   for  a  subregion   i,  S±      corresponds  to  H^,   C 
corresponds  to  Dlt      T   corresponds  to   (Dj  n  H^   and   F 
corresponds  to   (D]_  n  H0) .   For  generality,  suppose  true  contacts 
do  not  localize  a  target  to  a  single  subregion.   Then  after  a 
search  of  a  region  that  has  resulted  in  a  contact,  the  target's 
position  distribution  can  be  modified  as  follows: 

(27)  P(SjJC)  =  P(SjjT)  P(T|C)   +   P(SjjF)  P(F|C) 

where   i  =  1,2,  •••  ,n.   Eguation  27   results  since   T  =  T  n  C, 
F  =  F  n  C   and   P(SjJc)  =  {P[(Sj_  n  T)  +  P(Sj_  n  F)  ]  }/P(C)  .   The 
probability  p  =  P(t|c)   has  been  called  the  credibility  of  the 
contact.   In  terms  of  p  Eguation  27   becomes: 

(28)  P(Si|T)  =  P(Si|T)p  +  P(Si|F)(l  -  p) . 
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For  the  model,   P(S-jjT)   can  be  determined  by  the  coverage 
characteristics  of  the  detection  system  that  is  used  to  make  the 
contact.   In  particular,   P(SjjT)  =  P(T| S^) P(Si)/P(T)   where 
P(T)  =  2P(T| Sj ) P(Sj)   and  the  sum  index   j  =  1,2,  •••  ,n.   Then, 
using  the  correspondence  between  the  events   S^   and  H^   and  the 
events  T  and  (D±   n  H^   one  can  write   (Pd) i  =  p(Tlsi)   where 
(Pd) i  ^s  an  average  probability  of  detection  over  the  itn 
subregion.   In  one  positive  information  model,   PfSjjF)  =  P(Sj^)  , 
the  probability  the  target  is  in  the  itn  subregion  before  the 
contact  information  was  available,  and  p   is  determined 
subjectively  based  on  factors  associated  with  a  search.   The 
probability  determined  by   1  -  p  has  been  called  the  false  alarm 
probability.   However,  it  is  not  the  probability  of  a  false  alarm 
as  it  is  defined  in  Section  I.   Since   1  -  p  =  P(f|c),  it 
corresponds  to   P(H0|D1)   not  to   pf  =  P(D1|h0). 
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XV.   Search  Models  and  Search  Theory 

Search  theory  provides  a  basis  for  determining  optimal 
search  plans  for  a  target  whose  motion  and  location  are 
determined  within  some  bounds.   Here,  an  optimal  search  plan  is 
one  for  which  the  probability  of  finding  a  target  within  a  given 
length  of  time  is  a  maximum,  the  expected  time  to  find  a  target 
is  a  minimum  given  the  target  is  found  or  a  search  plan  for  which 
some  other  optimal  search  criterion  is  satisfied. 

Search  theory  results  are  based  on  models  of  the  search 
process.   To  the  degree  that  a  search  model  describes  a  search 
process,  an  optimal  search  plan  for  a  target  that  is  based  on  the 
search  model  should  provide  guidance  for  the  development  of  an 
operationally  feasible  search  plan.   However,  because  of  the 
limitations  of  analytical  search  models,  an  optimal  search  plan 
that  is  based  on  an  analytical  search  model  may  give  only  initial 
guidance  in  this  regard.   The  optimal  search  plans  that  are 
described  below  illustrate  this.   The  search  plans  are  based  on 
the  random  search  model.   Because  of  this,  the  requirement  on  the 
location  of  search  track  segments  is  not  realizable  and  the  time 
to  resolve  false  alarms  is  ignored. 

Optimal  search  plans  based  on  search  models  implemented 
through  a  monte  carlo  simulation  are  not  considered  here. 
However,  with  sufficient  information,  such  plans  have  the 
potential  of  being  both  implementable  and  more  optimal  in  a  real 
sense  than  an  optimal  search  plan  based  on  an  analytical  search 
model . 
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Three  Optimal  Search  Plans:   The  three  optimal  search  plans 
differ  through  their  definition  of  optimality.   However,  each  one 
is  based  on  the  following  search  model:   A  target  is  fixed  at 
some  point  in  a  region  that  consists  of   n   subregions.   A  search 
in  a  subregion  is  a  random  search  in  the  sense  of  the  definition 
in  Section  X   and  a  searchers  sweep  width  there  is  a  constant. 
In  addition,  a  search  of  a  subregion  will  not  detect  a  target 
which  is  in  another  subregion.   To  determine  a  plan,   let 
S^  =  (the  target  is  in  subregion  i}   for   i  =  1,2,  •••  ,n   and 
let  Pi  =  P(S^)   be  the  prior  probability  that  the  target  is  in 
the  i^-h  subregion.   Let  Wj_   be  the  sweep  width  in  the  i*-" 
subregion.   Let   6^  =  A^/W^   where  Aj_   is  the  area  of  the  i*-n 
subregion  and   5-j_   is   the  expected  track  length  to  find  the 
target  by  a  search  of  the  ith  subregion  given  the  target  is  in 
the  i*-*1  subregion,  a  characteristic  length.   The  probability   P 
that  the  target  will  be  detected  by  a  random  search  is  given  by: 
(29)       P  =  S  [1  -  expC-li/SiH-pi 

where   the  sum  index   i  =  1,2,  •••  ,  n   and   lj,   is  the  track 
length  of  the  search  in  the  i*-n  subregion. 

The  first  optimality  criterion  is:   Choose   1^   so  that   P 

is  a  maximum  subject  to  the  two  constraints:   1.  1  =  S  lj_   and 

2.  lj_  >  0   where  the  index   i  =  1,2,  •••  ,n.   Determining  this 

choice  is  a  nonlinear  optimization  problem  whose  solution  is 

given  in  Reference  19.   It  is: 

it/ Si   =  Infpi/Si)  -  L(k)    i  =  1,2,  •••  ,  k 
(30) 

1i/6i  =  °  i  =  k+l/k+2,  •*•  ,n 
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where   L(k)  =  (1/E6 j) «E[5 j «ln(pj/5 j) ]  +  1/ZSj   and  the  sum  index 
j  =  1,2,  •••  ,  k,  where  the  subregions  are  relabeled  so  that  the 
following  order  relation  holds:   P1/&1  >  P2/52  >  •••  >  Pn/^n  anc* 
where   k   is  chosen  so  that  for  k+l   the  solution  for   ljc+i 
using  L(k+1)   is  either  negative  or  zero. 

The  second  optimality  criterion  is:   Choose   lj^   so  that  P 
is  a  maximum  subject  to  the  two  constraints:   1.  c  =  2  c^   and 
2.  c^  >  0  where  the  index   i  =  1,2,  •••  ,n,  c^  =  kj/li   is  the 
cost  of  the  search  in  the  itn  subregion  and   k^   is  the  cost  per 
unit  track  length  in  that  subregion.    For  this  criterion,  the 
solution  to  the  corresponding  nonlinear  optimization  problem  can 
be  obtained  from  Equation  30  by  replacing   5^  by   e^  =  kj/S^ 
and  labeling  the  subregions  so  that   Pi/£i  >  P2/ei  >••••>  Pn/ei* 
The  basis  for  this  can  be  seen  by  replacing   li/&i   by   its 
equivalent   Cj/e^   in  the  exponential  term  in   Equation  29. 

The  third  optimality  criterion  is:   Choose   lj^   so  that  the 
expected  utility  of  finding  the  target  is  a  maximum  subject  to 
the  two  constraints:   1.   1  =  S  1^   and   2.   1^  >  0  where  the 
index  i  =  1,2,  •••  ,n.   For  this  criterion,  the  solution  to  the 
corresponding  nonlinear  optimization  problem  can  be  obtained  from 
Equation  30  by  replacing  p^  by  q-^  where  q^  =  Uj/p^   and  u^ 
is  the  utility  of  finding  the  target  given  it  is  in  the  ith 
subregion.   And,  in  addition,  labeling  the  subregions  so  that 
qi/&l  >  q2/&2  >  •••  >  <3n/5n'   Tne  basis  for  this  can  be  seen  by 
multiplying  the  summation  term  in  Equation  29   by   u^   so  that 
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Equation  29   gives  the  expected  utility  of  finding  the  target 
with  the  search. 

Equation  30   can  be  used  to  determine  an  order  of  search  for 
the  subregions  which  will  effectively  minimize  the  expected  track 
length  required  to  detect  a  target  given  it  is  detected.   To  do 
this,  divide  the  available  track  length   1   into  units  small 
enough  so  that  with  a  single  unit  only  the  1st  subregion  would  be 
searched.  Then  allocate  one  unit  to  the  search  of  the  Is*- 
subregion.   If  the  search  is  unsuccessful,  determine  the  optimum 
allocation  for  two  units.   Then  search  with  a  second  unit  so  that 
the  first  search  with  the  first  unit  plus  the  second  search  with 
the  second  unit  satisfy  the  optimum  allocation  for  two  units.   If 
the  search  is  unsuccessful,  continue  in  this  fashion  until  either 
the  target  is  found  or  all  the  track  length  is  expended.   That 
this  allocation  order  will  effectively  minimize  the  expected 
track  length  required  to  detect  a  target  given  it  is  detected  can 
be  argued  as  follows:   Let   L  be  the  track  length  at  detection, 
let   lu   be  a  unit  of  track  length  and  let   n  be  the  number  of 
units.   Then  the  value  of  the  probability   P(L  <  i*lu)   that  the 
target  will  be  detected  on  or  before  the  itn  step  of  the  search 
for  the  given  allocation  order  will  be  greater  than  or  equal  to 
its  value  for  any  other  allocation  order  with  the  same  allocation 
step  size.   Since  the  value  of   P(L  <  1)   will  be  equal  to  its 
value  for  any  other  allocation  order  of  the  optimum  allocation 
and  since   P(L  <  i*lu|L  <  1)  =  P(L  <  i»lu)/P(L  <  1),  the  value  of 
the  distribution  function  FL(i»lu|L  <  1)  =  P(L  <  i*lu|L  <  1) 
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will  also  be  greater  than  or  equal  to  its  value  for  any  other 
allocation  order.   This  implies  that  the  expected  track  length 
given  detection   E(l|l  <  1)  =  2  [1  -  FL(i»lu|L  <  1)]   where  the 
sum  index   i  =  1,2,  •••  ,  n   is  effectively  a  minimum  for  the 
given  allocation  order.   A  search  based  on  the  optimum  allocation 
given  by  Equation  30   and  the  given  allocation  order  is 
equivalent  to  the  following  search  prcedure:   After  an  allocation 
of  track  length   lu   and  an  unsuccessful  search,  new  values  for 
P(S^)   are  calculated  using  Equation  26   and  then  Equation  30   is 
used  with  these  new  values  to  determine  the  next  optimum 
allocation.   A  discussion  of  this  procedure  is  given  in  Reference 
6.   And  an  example  of  its  application  is  given  in  Reference  20. 

Equation  30   also  defines  an  optimal  search  plan  for  a 
detection  system  that  searches  beams  and  can  be  described  by 
Equation  28   by  replacing   lj_   by  t^  where   t^   is  the  time  the 
itn  beam  is  searched  and  by  replacing  8^      by   t^   where   tj_   is 
the  expected  time  to  detect  the  target  by  a  search  of  the   ith 
beam  given  the  target  is  in  the  i*1*1  beam,  a  characteristic  time. 

For  a  more  extensive  discussion  of  search  theory  and  its 
application  to  military  operations  research,  see  Reference  21. 
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